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ABSTRACT 

We present the observations of GRB090510 performed by the Fermi Gamma- Ray Space Telescope 
and the Swift observatory. This is a bright, short burst that shows an extended emission detected 
in the GeV range. Furthermore, its optical emission initially rises, a feature so far observed only in 
long bursts, while the X-ray flux shows an initial shallow decrease, followed by a steeper decay. This 
exceptional behavior enables us to investigate the physical properties of the GRB outflow, poorly 
known in short bursts. We discuss internal shock and external shock models for the broadband energy 
emission of this object. 
Subject headings: gamma rays: bursts 



1. INTRODUCTION 

With the availability of a relatively large sam- 
ple of Gamma-Ray Bursts (GRBs), we came to 
recognize that they com prise of two large classes 
(|Kouveliotou et al. 19931 ): the so-called short-hard 
GRBs (duration < 2 s) and the long-soft ones (> 
2 s). There is now increasing consensus that the ob- 
served dichotomy among long/short GRBs may indi- 
cate diverse initial physical conditions and progenitors. 
Long GRB s are associated w ith the demise of mas- 
sive stars (jFerrero et al. 20061) . Ins tead, short GRBs 
often occur in ear l y-type galaxies ( Zhang et al. 20091 : 
iGehrels et al. 2009t iFong et al. 20091 1 This supports 



At 00:23:00 UT, May 10th, 2009, BAT, which oper- 
ates in the 15-350 keV range, triggered on GRB090510 
(|Hoversten et al. 2009f ) . Swift slewed immediately to the 
burst. The duration was T 90 = 0.30 ± 0.07 s. Detailed 
analysis of the BAT data is shown in Ukwatta et al. 
(2009). 



2.2. XRT data 



The 



Swift 



X-ray 



their interpretation in terms of compact object mergers. 

Crucial information on GRBs is revealed by 
their afterglows, w hich can be monitored by Swift 
(|Gehrels et al. 20041 ) in the optical and the X-ray range 
as soon as ~ 100 s after the burst. In this paper, we 
present the study of the short GRB090510 with Swift 
and Fermi in a broad energy range, which extends from 
the optical up to a few GeV. 

We report our observations and analysis in §2; in §3 
we propose two different interpretations, in §4 we draw 
our conclusions. Hereafter, we use the conventions X = 
10 n X n for cgs units and F cx t~ a v~ 13 , where F is energy 
flux, t is t he time from the trigger of Sw ift Burst Alert 
Telescope (|BAT: Barthelmv et al. 2005af) and v the fre- 
quency. Errors are reported at la, unless otherwise spec- 
ified. We assume a cosmology in which Ho = 71 km s _1 
Mpc" 1 , Q m = 0.27, Q x = 0.73. All the fluxes, times and 
frequencies are measured in the observer's frame. 

2. OBSERVATIONS AND DATA ANALYSIS. 
2.1. BAT data 



* Corresponding authors: M. De Pasquale: mdp@mssl.ucl.ac.uk; 
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Telescope 

(|XRT: 0.3-10 keV: Burrows et al. 2005T) began to 
observe the X-ray afterglow of GRB090510 at T+98 s. 
The lightcurve (Fig. 1) shows an initial slow flux decline. 
Observations were interrupted when the source entered 
the Earth constraint at T+1.9 ks. When they resumed, 
at T+5.1 ks, the flux was much lower. A broken pow- 
erlaw fit of the lightcurve gives as best fit parameters 
an early decay slope ax,i — 0.74 ± 0.03, break time 
t x = 1.43±g;5§ ks, late decay slope a x ,2 = 2.18 ± 0.10; 
X 2 = 112 with 77 degrees of fredom (d.o.f.), which is still 
marginally acceptable (chance probability P = 0.0054). 



2.3. UVOT and other optical data 



The 



Swift Ultra Violet and Optical Tele- 
scope (IUVOT: 160-800 nm: Roming et al. 20051 : 
iPoole et al. 2008t ) began settled exposures at T+97 s. 
The lightcurve of the optical afterglow, produced by 
renormalizing all individual filters to white, as described 
in Oates et al. (2009), is shown in Fig. 1. The op- 
tical emission rises until ~ 1.6 ks, then decays. The 
optical lightcurve is well fitte d (y 2 = 23.9 with 19 
d.o.f.) by a broken powerlaw (jBeuermann et al. 19991) 
with a smooth break. The best fit parameters 
0-Opt,i = — 0-50i 013 ; t pea k 



are: 



aopt.2 



1.13 



+0.11 

-o.io- 



— 1.58^0.37 ks; 



Adding a constant does not 



improve the fit significantly, suggesting a small host 
galaxy contribut ion. Very Large Telescope observations 
(|Rau et al. 2009?) provide a spectroscopic redshift of 
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z = 0.903. Using this redshift and the Fermi spectral pa- 
rameters, the isotropic equivalent energy of GRB090510 
is E iso = 1.08 x 10 53 erg s in the 10 keV - 30 GeV rest 
frame (|Abdo et al. 2009f ). 

2.4. Fermi data 

GRB090510 triggered both instruments on- 
board the Fermi o bservatory (jGuiriec et al. 20091 : 
lOhno. fc Pelassa 2 009). The Gamma-ray Burst Moni- 
tor (GBM; 8 keV - 40 MeV) observed the burst during 
the prompt emision phase, and after an autonomous re- 
pointing, the Large Area Telescope (LAT; 20 MeV - more 
than 300 GeV) began observations and detected a long- 
lasting (up to 200s) high-energy (up to 4 GeV) emission. 
The analysis and interpretation of the prompt emission 
is presented in Abdo et al. (2009). Follow up observa- 
tions lasted until 1500 s when the source was occulted by 
the Earth, and resumed ~ 3.5 ks later. The observation 
epochs are defined in Table [TJ A time-resolved spe ctral 
analysis using diffuse events (jAtwood et al. 20091 ) was 
performed (Table [T]). The spectrum shows no significant 
evolution, and it is well-fitted by a powerlaw with 
energy index /3 7 = 1.1 ± 0.1. The lightcurve shows no 
significant feature and is best fitted by a powerlaw. The 
onset of the GBM emission (T+0.013 s) is a sensible 
reference for this temporal fit and yields a decay index 
OLy = 1.38 ± 0.07 (x 2 /dof = 9.4/7) (Fig. 1). For the 
binned spectral analysis described in §3.2, transient 
events from an energy-dependent region of interest were 
used, and front and back events were treated separately 
(|Atwood et al. 2009D . 

3. DISCUSSION 

GRB090510 was a short burst with a relatively bright 
afterglow, and Ej SO is among the highest for this class 
(jGraham et al. 20091 ). The early rise of the optical flux 
is so far unique in short GRBs. More importantly, 
GRB090510 shows high energy emission up to the GeV 
range, until T+200 s. 

Energetic s hort GRBs with optical transients, such as 
GRB050724 (jBarthelmv et al. 2005bl) , typically have ex- 
tended emission (EE) detected b y BAT and XRT fol- 
lowing the hard emission spike jTroia et al. 20081 ). If 
GRB090510 had occurred at z = 0.26, as GRB050724, it 
would have produced a flux in the BAT range of a few 
10 -9 ergs cm~ 2 s _1 until ~ 100 s, and we would have 
classified it as an EE-GRB. 

The nature of this high energy emission is nevertheless 
not easy to understand. EE often fades slowly for a few 
hundreds of seconds, then vanishes with a slope which 
can be as fast as a ~ 7; after this sudden drop a late af- 
terglow with a typical decay slope a ~ 1.4 is sometimes 
observed (e.g. GRB050724; GRB080123, Manganoet al., 
2008). The fast decay and the extrapolation of the late 
afterglow back to early epoc hs suggest th at EE is not the 
onset of the late afterglow (|Nakar 20071 ). Furthermore, 
in a few case s where the EE is bright enough to be stud- 
ied in detail (jNorris fc Bonn ell 2006) , it is found to have 
variatio ns too rapid to be expla ined with external shock 
models (Mesz aros fc Rees 19 93). EE might instead in- 
dicate a declin i ng activity of the GRB central engine 
(iRosswog 20071: iMetzger et al. 20081: iPerna et al. 20061 : 
IGoad et al. 2007] ): once this activity ends, then falls 
abruptly, and the forward shock (FS) emission prevails. 



In other cases, however, the flux decay from the begin- 
ning of Swift observations seems due to the usual FS 
mechanism, such a s in GRB051221 flB urrows et al. 20061) 
and GRB 061201 (|Stratta et al. 20071) . 

We propose and discuss separately 2 scenarios to ex- 
plain the emission after the initial spike: in the first one, 
the emi ssion is due to both exte rnal FS and internal 
shock, (|IS. Rees fc Meszaros 1994|) while in the second 
the emission is due to FS alone. 

3.1. X-ray internal shock, optical external shock 

In the first scenario we assume that the initial X- 
ray (until the break at ~ 1.4 ks) and 7-ray flux are 
IS emission, while the FS is responsible for the opti- 
cal lightcurve and the late X-ray flux. In particular, 
the optical rise may be due to the onset of FS emis- 
sion, detected 10 2 — 10 3 s after the trigger in long GRBs 
(jOates et al. 20091 ). This model can explain the different 
behavior of the early X-ray/LAT and optical lightcurves. 

We constrain some physical properties of the FS 
blastwave. The initial Lorentz factor of the ejecta is 

T = 1.4 x 10 2 K^ 8 n-^H'^lz (|Sari 19971 ) . where E 
is the isotropic kinetic energy, n the environment den- 
sity in cm -3 , and t pea k is the peak time. The max- 
imum FS flux is at the synchrotron characteristic fre- 
quency u m , and it is F Vm = 1.3 x lQ^E^e 1 ^ 2 _ 2 ?i 1/ ' 2 /iJy 
(|Granot fc Sari 2002f) . where e# is the fraction of inter- 
nal energy in the magnetic field. 

Th ese parameter valu es must be consistent with Tq > 
1000 (jAbdo et al. 20091 ) and with the UVOT data, which 
give a 3cr lower limit on t pea k > 730 s and a peak flux F ~ 
100/iJy. The first constraint can be written E^n^ 1 > 
2.6 x 10 6 . If v m is just below the optical band, the con- 
straint on the flux becomes E^e^^n 1 / 2 ~ 7.7 x 10 -3 . 
Assuming £3,-2 — 1, the model is consistent with obser- 
vations for E 53 > 5.4, while n ~ 5.9 x lO -5 ^ 2 . 

The XRT and LAT fluxes can be explained by IS 
synchrotron emission with v c < vopt < v sa < vx < 
v m , where v c is the synchrotron cooling frequency 
and v Ra is the synchro tron self-absorption frequency 



(|Guetta fc Granot 2003[ ). The synchrotron luminosity 
estimated from the 100s spectral energy distribution 
(SED, Fig. 2), is L ~ 10 50 erg s" 1 . We find that for 
this value of L and for p = 2.4, e ej _i = 5.5, £3,-2 = 33, 
T = 410, t v = 3 x 10~ 5 s, where p, e e , T and t v 
are the index of the powerlaw electron energy distri- 
bution, the fraction of energy given to electrons, the 
bulk Lorentz factor and the variability timescale respec- 
tively, we have v m ~ 210 keV, F(1.7 keV) ~ 75 //Jy, 
which is within a factor ~ 3 from the observations, and 
F(100 MeV) ~ 4.2 x 10~ 3 fiJy, which is consistent within 
2cr of the data. The cut-off energy for pair production 



is hi 



77 



1.6 GeV, thus allowing the late emission of 



~ 1 GeV photons. IS does not produce detectable emis- 
sion in the optical, since this is below v sa ~ 0.32 keV. 

Compared with the scenario presented in §3.2, this 
model has the advantage of not requiring extreme val- 
ues of Tq to explain an early (few seconds) FS emission 
onset in a low density environment expected for a short 
burst. However, it needs some fine tuning of parame- 
ters. The optical rise slope ao P t.i — —0.5 is shallower 
than that expected at the FS onset (a — —2), although 
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similar slow rises have been observed ([Oates et al. 20091) . 
It is possible that the onset of our observations caught 
the end of this steep rise phase, when the afterglow was 
turning over to a decay. Another possible problem is 
that the required density appears very low. We note 
that the observed slope would be expected if v m were 
crossing the optical band and, in general, broad FS on- 
set rises are also expected for outflows observed off- axis 
(|Panaitescu &i Verstrand 20081 ). However, a bright and 
hard event such as GRB090510 is difficult to reconcile 
with the latter s cenario, which predic ts soft and dim 
prompt emission (jYamazaki et al. 2002D . 

3.2. Optical, X-ray and GeV emission from external 

shock. 

A second possibility is that the afterglow of 
GRB090510, including the emission detected by LAT, 
is entirely produce d by the FS prop agating in a constant 
density medium (|Sari et al. 19981 ). According to the 
model, the broad afterglow spectrum 1 consists of three 
segments: a low-energy tail, of spectral slope fa = —1/3; 
another segment, for v m < v < v c , where fa = (p — l)/2; 
blueward of v c , the third segment has fa = p/2. For 
comparison with this spectral template, we built up 5 
SEDs, at 100 s, 150 s, 1 ks, 7 ks, and 12 ks (Fig. 3), 
all including UVOT and XRT data, and LAT data were 
also included in the first SED. LAT data were accumu- 
lated between 10 s and 200 s (i.e. well after the end 
of the prompt emission seen in the GBM), and renor- 
malized to 100 s using the decay index of a 7 = 1.38. 
We fitted the SEDs simultaneously with a double broken 
powerlaw model, forcing fa = —1/3 and fa = fa + 1/2. 
We allowed the breaks to vary and Galactic and host 
extinction were accounted for. The result is acceptable 
(x 2 /d-o.f.=110.3/83) and is shown in Fig. 3 and Tabled 
The FS alone could successfully describe the spectrum 
over 9 decades of frequency. A break between X-ray and 
7-ray ranges is fitted, at E\ ~ 300 MeV, but not con- 
strained. It is studied more precisely by fitting the 100 s 
SED alone, freezing N H and E(B - V) at the 5 SEDs fit 
results and leaving fa and fa free to vary (see Tabled- 
This fit fulfills the relation fa = fa + 1/2 (at 1.3 a) and 
a significant break is found (3.6 a), although it yields a 
slightly harder fa (1.8 a) than the 5 SED fit shown in 
Fig. 3. The LAT emission shows no spectral evolution, 
even at early times. Therefore, to better characterize the 
high-energy spectrum, the SED at 100 s was rebuilt in- 
cluding LAT data between 0.38 s and 200 s (see Tabled] 
and Fig. 2). A significant break (> 4.5c including sys- 
tematics) between 10 and 133 MeV was found. However, 
including this selection of LAT data in the 5 SED fit 
yields a worse fit (x 2 /dof = 125.3/83) than that shown 
in Fig. 3. 

In this FS interpretation, the initial increase of the op- 
tical emission is due to v m approaching the optical band. 
The X-ray is already decaying because it lies above v m . 
In order to verify whether the required physical parame- 
ters are plausible, we impose the following constraints: i) 
F(l ks) ~ 2.2 fiJy at 10 18 Hz and ii) v m {\ ks) ~ 10 16 Hz. 
Adopting the expressions for F Um , v m , and v c from 
Granot & Sari (2002), the constraints above lead to 

1 The self-absorption frequency is not relevant in this study. 



the following equations: e_B,-2 — 14 E 5 ^' 2 e e 4 _ 1 ^ p 4 'Vg S , 
n ~ 1.5 X 10~ 6 (1330) p - 2 5 £; 53 1 e4_ i ^4 I/ - j p-i j where 

= 3(p - 2)/(p - 1) and v os = u m (l ks)/10 16 Hz. We 
verified that the synchrotron self-Compton cooling is not 
significant for t < 1.5 ks and v c ( 1.5 ks) > 10 18 Hz for 
a reasonable range of parameters (|Nakar et al. 20091 ). A 
very low, but not implausible, density is suggested. 
The flux at lOOMeV is 

F v>Uc ~ 2.4 x 10- 3 (2.2 x 10- 3 )(P- 2 - 5 )£ 53 e e ,-i 

(i/100s)( 2 - 3 P)/ 4 (/i^/100MeV)-P/ 2 ^^r 2)/ Vjy 

(1) 

This is consistent with the LAT data at 100s, provided 
E^te,-! — 5, p w 2.5, £ p w 1, and v os w 1. For these 
parameters, v c <C 4GeV at t > 1 s, so that the flux in the 
LAT energy range is approximated to be F oc £( 2_3 p)/ 4 ~ 
t~ 1A at t > 1 s, consistent with the LAT light curve. We 

note that r > dSOOE^nZ^ 8 is required for the FS 
onset time to be < 1 s. 

In summary, the spectral properties of GRB 090510 
could be explained by a simple FS model. We note, 
though, that this simple model is hard to reconcile 
with some of the observed temporal properties. Firstly, 
it predicts an X-ray decay index before the break of 
a = 3 fa/2 = 1.16 ± 0.06, clearly inconsistent with the 
observed ax : i = 0.74±0.03. Secondly, i f the X-ray break 
at t = 1.4 ks is attributed to a jet break (|Sari et al. 19991) 
and, after this time, optical and X-ray lie on the same 
spectral segment, then the asymptotic optical decay in- 
dex should be the same as ax, 2 = 2.18 ± 0.10. How- 
ever, a fit of the whole lightcurve with a smooth broken 
power-law (§2.3), adding a constant (as a host galaxy 
contribution), gives an asymptotic decay slope aopt.2 = 
1.131 

at)9' incompatible with the X-ray decay. Finally, 
although the error bars are quite large, we notice that, 
taken at face value, the slope of the UVOT spectrum in 
the 1 ks SED is negative (Fig. 3), suggesting that v m 
may be already below the optical at that epoch. 

The above mentioned flaws imply that the simple 
FS model is not viable to explain the properties of 
GRB090510. However, this model relies on highly ide- 
alized assumptions, and it is known that several GRB 
afterglows do not strictly follow its simple predictions. 
Plausible effects that may affect the predictions and ease 
the comparison with GBR09051 are: 

- a phase of energy injection (jSari fc Meszaros 2 000). 
or an evolu tion of the microphysic al parameters of the 
blast wave ([Panaitescu et al. 2006f) : both may cause an 
early shallow decay of the X-ray flux; 

- the transit of v m slightly after the jet break, which 
could explain a shallow late optical decay. 

With the present data, we are unable to distinguish 
between energy injection and microphysical parameter 
evolution. As for the X-ray decay post jet break, hydro- 
dynamical simulations show that in itially the jet decay 
slope can reach a ~ 3 ()Granot 20071) . Both the processes 
of energy injection and parameter evolution are capable 
of keeping the decay shallower, so that a late X-ray decay 
slope of ax, 2 = 2.18±0.10 could be achieved. Therefore, 
with some extensions, the FS model could arrange the 
temporal properties, although some fine tunings would 
be needed. 
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4. CONCLUSIONS 

We have reported the Swift and Fermi observations 
of the short GRB090510, an event endowed with bright 
prompt and afterglow emission, and detected in the GeV 
range up to 200 s after the trigger. The initial X-ray 
emission shows a slow decay up to ~ 1.4 ks, after which 
it quickly drops. The optical flux peaks at ~ 1.6 ks. We 
have explored two scenarios to explain the observed be- 
haviors. 

In the first scenario, the early flux detected by XRT and 
LAT is due to IS, while the optical rise is the onset of 
FS emission or the transit of v m . This interpretation 
does not require extremely high values of Lorentz factor, 
should the density of the environment be very low. We 
also find that reasonable values for the physical param- 
eters can lead to the observed properties, which might 
favour the model, although some fine tuning is neces- 
sary. The second scenario assumes that the FS produces 
the full spectrum of the emission, observed from the op- 
tical to the GeV band. The 7-ray, X-ray and optical 
spectrum can be reproduced by the template FS spectral 
model and the required physical parameters are plausi- 
ble. Although the simple FS model fails to reproduce 
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°- 85 +0.30 
1 20" ' 20 
fl qcj-0.26 

i .,-,-0.28 
L41 +0.31 


10- 2 

imtltt 10- 2 

c 7+2.4 ln _3 
O. i _ j g ±U 

b - 4 -1.6 iu 




5: 2.5 - 5.5 
6: 5.5 - 11.5 
7: 11.5 - 37. 
8: 37. - 69.5 
9: 69.5 - 200. 


163. 
58. 
71. 
59.9 
43. 


°- 76 +0.26 

0.86 +0 ; 44 

o 97— 0.59 
z - z '+0.70 

u -°°+0.39 

L - ' 4 +0.71 


a +3.6 -i n — 4 

8-4_2 7 10 
2.0+Q 9 10_4 
1.67« ; j 2 10- 4 
AA+_\l IO" 5 
i.6±^;° 10" 5 




10: 200. - 400. 


~ 0. 


1.1 / 0.5 / 2.5 


< 4.7 / 3.2 / 7.9 10" 


6 


11: 400. - 800. 


~ 0. 


1.1 / 0.5 / 2.5 


< 2.3 / 1.6 / 3.9 10" 


6 


12: 800. - 1500. 


~ 0. 


1.1 / 0.5 / 2.5 


< 2.1 / 1.2 / 3.8 10" 


6 


13: 4200. - 7200. 


~ 0. 


1.1 / 0.5 / 2.5 


< 0.48 / 0.3 / 1.0 10 


-6 


14: 10150. - 13000. 


~ 0. 


1.1 / 0.5 / 2.5 


< 0.47 / 0.3 / 0.9 10 


-6 


15: 15800. - 18500. 


~ 0. 


1.1 / 0.5 / 2.5 


< 0.52 / 0.3 / 1.0 10 


-6 



TABLE 1 

lat time-resolved spectroscopy. the source test statistic was defined as twice the difference of the unbinned 
log-likelihood between the null hypothesis (background only) and the alternative hypothesis (presence of a source) . 
Bottom table shows 95% confidence level upper limits on flux for different energy indeces. 



SED 




5 SED 


100 s 


100 s 


LAT dataset (s) 




10. - 200. 


10. - 200. 


0.38 - 200. 


N H (x 10 21 cm 2 ) 




1 ^9 +0.33 
± - oz -0.30 


1.52 (fixed) 


1.52 (fixed) 


E(B - V) (mag) 




000 + - 005 
u.uuu _ .ooo 

r 0.43l»' ° (100 s) 


0. (fixed) 


0. (fixed) 






0.17l»; ( « (150 s) 


01 +0.05 
-0.06 


n 01 +0.06 


El (keV) 


i 


0.0371°;"? (1000 s) 
< 0.001 (7000 s) 
. < 0.01 (12000 s) 


ft 




0.77 ± 0.04 


61 + 06 


62 + ' 08 


E\ (MeV) 




~ 300 (100 s) 


[20 - 135] 


[10 - 133] 


03 




ft + 1/2 


1 44 +0-26 
J-- 44 -0.22 


1 14 +0-10 
J-- 14 -0.09 



TABLE 2 

Best fit parameters obtained by fitting the 5 SEDs simultaneously or the SED at 100 s only with a double broken 
powerlaw model (see text). Nh and E(B — V) ARE host absorption and extinction RESPECTIVELY; E\ and E\ are the two 

break energies calculated at the epoch in parenthesis. 



G 



De Pasquale et al. 




Time since BAT trigger (s) 



Fig. 1. — Top: LAT flux above 100 MeV and best fit to the flux decay (line). Bottom: energy flux densities averaged in the observed 
energy bands: BAT (15 keV - 350 keV, stars); XRT (0.2 keV - 10 keV, crosses); UVOT renormalised to white (diamonds); LAT (100 
MeV — 4 GeV, filled squares; the average spectral index was used to convert from photon to energy flux) with upper limits for /3 = 1.1 
(triangles). The prompt emission is shown for comparison: GBM (8 keV — 1 MeV, circles), LAT (100 MeV - 4 GeV, empty squares). XRT 
lightcurve is obtained as in Evans et al. (2007, 2009). 

All data are shown with 68% error bars or 95% confidence level upper limits. 



Early afterglow of the short GRB090510 
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Fig. 2.— UVOT- 



-XRT-LAT count spectrum, with best fit and residuals shown (see text) 
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Fig. 3.— UVOT-XRT-LAT SEDs at different epochs, with the best fit shown (see text). The butterfly at 100s indicates the 68% 
confidence level region for the LAT flux, obtained from an unbinned likelihood analysis (95% error bar at 100 MeV is shown). Successive 
SEDs in time order are rescaled by 1:1, 1:10, 1:100, 1:1000, 1:10000. 



